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Suspensions of motile cells are model systems for understanding the unique mechanical properties of living materials which often 
consist of ensembles of self-propelled particles. We present here a quantitative comparison of theory against experiment for the 
rheology of such suspensions. The influence of motility on viscosities of cell suspensions is studied using a novel acoustically- 
driven microfluidic capillary-breakup extensional rheometer. Motility increases the extensional viscosity of suspensions of algal 
pullers, but decreases it in the case of bacterial or sperm pushers. A recent model [Saintillan, Phys. Rev. E, 2010, 81, 56307] 
for dilute active suspensions is extended to obtain predictions for higher concentrations, after independently obtaining param¬ 
eters such as swimming speeds and diffusivities. We show that details of body and flagellar shape can significantly determine 
macroscale rheological behaviour. 


1 Introduction 

Many living materials, such as suspensions of motile microbes 
or of ATP-powered cytoskeletal polymers, are large ensembles 
of nearly identical and motile subunits that interact strongly 
with their neighbours. Understanding the properties of such 
systems presents unique conceptual challenges. Since each 
elemental subunit such as a motile cell is by itself a driven- 
dissipative system, their collectives operate well out of equi¬ 
librium even in the absence of any external forcing. Local 
fiuctuations in motion further have a non-thermal origin, and 
the fiuctuation-dissipation theorem has been shown to be in¬ 
applicable in these systems .1^ 

The recent development of a continuum framework that 
shows that these intrinsically non-equilibrium systems share 
universal features is therefore a significant theoretical break- 
through.l^ This theory of “active matter” suggests that a net 
average local alignment of self-propelled particles must, in a 
continuum description, lead to a contribution to the stress ten¬ 
sor arising from propulsive forces or “activity” of the particles. 
An interesting prediction is that particle activity must change 
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the macroscopic viscosity of the suspension. An axisymmetric 
self-propelled particle exerting a net propulsive thrust along 
its principal axis generates a hydrodynamic force dipole. The 
flow field around a single E. coli cell has been measured to be 
approximately that of a positive hydrodynamic dipole® which 
forces the ambient fiuid around each particle axially outward 
towards its two ends. The suspension viscosity for such “push¬ 
ers” is predicted to decrease below its value for a passive sus¬ 
pension of inactive particles of the same size, shape and con¬ 
centration. Conversely, a suspension of “pullers” with nega¬ 
tive hydrodynamic dipoles will have a higher viscosity than 
a passive suspension. These predictions® are su ppo rted by 
experimental observations with bacterial pushers®^ and al¬ 
gal pullers.® The observations appear to clearly confirm a key 
generic feature of the rheology of active suspensions. 

Quantitative micro structural models relating particle size, 
shape, concentration and motility to rheological properties are 
just beginning to emerge,®^ but thus far there have been no 
systematic comparison of their predictions against experimen¬ 
tal data. We present here such a comparison for suspensions 
of wild-type strains of the microalga Dunaliella tertiolecta, 
the bacterium Escherichia coli and mouse spermatozoa. D. 
tertiolecta uses its pair of flagella in a manner similar to the 
puller Chlamydomonas reinhardtii studied by Rafai et al.^ 
whereas E. coli and the sperm use flagellar tails to push for¬ 
ward, like B. subtilis studied by Sokolov and Aranson.®^'. coli 
cells have multiple flagella and use run-and-tumble swimming 
by bundling and unbundling their fiagella. Sperm were cul¬ 
tured under conditions to induce capacitation and hyperactive 
swimmingTrajectories of cells swimming were analyzed to 
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characterize their motility under quiescent conditions. Mea¬ 
surements of the extensional viscosity of live and dead cell 
suspensions were obtained with a novel microfluidic rheome¬ 
ter developed by.^ A recent model for the rheology o f di- 
lute suspensions of self-propelled hydrodynamic dipoles^®^ 
is extended here to the non-dilute regime to relate the non- 
Newtonian elongational viscosity of active suspensions to par¬ 
ticle volume fraction. We demonstrate here that despite the 
wide range of particle size, shape and motility characteristics, 
bulk mechanical behaviour of active suspensions may be accu¬ 
rately characterized in terms of a small number of parameters. 

2 Experiments 

2.1 Cell culture and suspensions 

D. tertiolecta Butcher was collected and isolated from Port 
Phillip Bay, Victoria, in December 2009. Cultures were 
maintained at 20°C in modifled F-medium (30 g/L aquar¬ 
ium salt, 250 mg/L NaNOs, 18 mg/L KH 2 PO 4 , 9 mg/L 
iron(III) citrate (CeHsOyFe), 9 mg/L citric acid (CeHgOy), 
0.2 mg/L MnCl 2 . 4 H 20 , 0.023 mg/L ZnS 04 . 7 H 20 , 0.011 
mg/L C 0 CI 2 . 6 H 2 O, 0.005 mg/L CUSO 4 . 5 H 2 O, 0.008 mg/L 
Na 2 Mo 04 . 2 H 20 , 0.65 ji g/L H 2 Se 03 and traces of vitamin 
B12, biotin and thiaminel^. 20 mL of inoculum were added 
to 400 mL of F-medium and incubated in 2-L glass bottles in 
a laboratory growth cabinet at 20 oC 0.1 °C with a 16:8 light 
dark cycle using white fluorescent lights with a photon flux 
of 60 /imol photons/ (m^ s). Cultures were bubbled with air 
through an aquarium air stone to provide a source of inorganic 
carbon (CO 2 ). After 5 days, F-medium was added to bring 
the total culture volume to 1.5 L. Samples for experimenta¬ 
tion were harvested after a further 6 days during the log-phase 
of growth into 50 mL polypropylene capped-tubes and cen¬ 
trifuged at 3500 rpm at 20 °C for 10 minutes to collect cell 
pelletl^. 

Wild-type E. coli K 12 strain was procured from ATCC, 
USA (# 10798). Standard media in the form of Luria Bertini 
(LB) broth (# L3022, Sigma Aldrich; 10 g/L tryptone, 5 g/L 
yeast extract, 5 g/L NaCl) and/or Luria Agar (# L2897, Sigma 
Aldrich) was used for bacterial culture. A UV-VIS spec¬ 
trophotometer (# UV-2450, Shimadzu) was used to charac¬ 
terize the bacterial growth by measuring absorbance / optical 
density at 600nm. About 0.5 mL of sterile LB broth was put 
into the sterile ATCC vial containing lyophilized culture. A 
small amount of the suspended culture (around 0.05 mL) was 
inoculated on to sterile Luria agar slants, and incubated at 37 
°C for 18-20 hours. A single colony from an agar slant was 
transferred to 5 mL of sterile LB broth and incubated at 37 
°C for 6-7 hours with vigorous shaking (at 170 rpm), till ab¬ 
sorbance at 600 nm reached 0.4. About 0.3 mL of 50 wt% 
glycerol-water solution was added to 0.7 mL of this mid-log 


phase culture and stored at -73 °C for future use. From the 
glycerol-freeze stock, a tiny amount is scraped off and inocu¬ 
lated under aseptic conditions to 5 mL of sterile LB medium. 
The culture was incubated at 37 °C for 16-18 hours with vig¬ 
orous shaking (at 170 rpm). A small amount (around 0.05 mL) 
was transferred into sterile LB media of 160 mL volume in a 
shake-flask. Cultures were incubated at 37 °C for 6-7 hours 
with vigorous shaking (at 170 rpm). 

Sperm from C57BL (wild-type) mice, extracted from cauda 
epididymii using the back-flushing method^ was added to 5 
mL pre-warmed MT 6 medium (125 mM NaCl, 2.7 mM KCl, 
1 mM MgCl 2 . 6 H 20 , 0.35 mM NaH 2 P 04 . 2 H 20 , 5.5 mM glu¬ 
cose, 25 mM NaHCOs, 1.7 mM CaCl 2 . 2 H 20 , 60 jiM bovine- 
serum albumin) containing methylcellulose and incubated at 
37°C for 90 min. 

Algal cultures were centrifuged at 3500 rpm for 8 min. 
at 20° C and pellets were re-suspended in growth medium at 
desired concentrations. LugoTs iodine (100 g/L KI, 50 g/L 
iodine crystals) was added to Eppendorf tubes to kill algal 
cells. These were centrifuged, the supernatant removed and 
re-suspended in growth medium again. E. coli cultures after 
log-phase growth were centrifuged at 4550 rpm for 10 min. 
at 4°C, and pellets were re-suspended after weighing in a (pH 
8.2) buffer of 10 mM K2HPO4, 0.1 mM EDTA and 0.2 wt.% 
glucose to prepare suspensions of various cell volume frac¬ 
tions. Suspensions were exposed to UV light for 30-60 mins, 
to kill cells without signiflcant lysis. Fresh sperm samples 
prepared as above were first tested in capillary-thinning exper¬ 
iments. Standing sperm suspensions for 90 min. inactivated 
motility. 


2.2 Particle tracking 

Image-analysis was used to determine the average fractional 
area covered by cells at the focal plane; this was assumed to 
be equal to the volume fraction. In bacteria and algae, flag¬ 
ellar filaments were not resolved in the image analysis. Sam¬ 
ple droplets were placed on teflon-coated glass slides, with 
a cover slip on top. The gap width between slide and cover- 
slip was typically 1 mm. Microscope videos for bacteria and 
algae were captured with 20X and lOOX (Olympus) lenses, 
respectively, and a high-speed camera (SA5, Photron; 50 fps; 
768 X 816 pixels). Swimming speeds and diffusivities were 
determined by processing images (Image!) and cell tracking 
(Imaris). Sperm suspensions were loaded onto both cham¬ 
bers on a 2X-Cel 80 jim slides chambers covered with 2X- 
Cel Cover Glass (Hamilton Thome Research). Slides were 
inserted into a Hamilton Throne IVOS for computer-aided se¬ 
men analysis (CASA). At least 1000 sperm were counted in 
each chamber. Sperm motility characteristics were analyzed 
through image analysis as done for the algae and bacteria. 
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2.3 Acoustically-Driven microfluidic extensional rheom- 
etry 

Unlike shear rheometry, techniques for reliable measurement 
of fluid properties in extensional flows have been established 
only relatively recently.^ In capillary-breakup extensional 
rheometry (CaBER), a liquid bridge is usually first created by 
rapidly moving apart end-plates between which a sample drop 
is sandwiched. If the end-plate separation is large enough, 
the bridge subsequently begins to thin due to the Rayleigh- 
Plateau instability. The rate at which a liquid bridge thins is 
governed largely by the balance of the capillary stress against 
the inertial and viscous stresses induced by the extensional 
flow about the necking plane and hence it is in principle pos¬ 
sible to extract the viscosity by monitoring the neck radius 
as a function of time.^ This technique has in the past been 
used for highly viscous samples. ^3 Obtaining reliable mea¬ 
surements with low-viscosity complex fluids such as aqueous 
cell suspensions however presents two challenges. Firstly, me¬ 
chanical motion of end-plates sets off inertial instabilities that 
quickly break up liquid bridges.^ The motion within the liq¬ 
uid bridge following the sudden stopping of the end-plates is 
complex and is not described by a simple stress balance. Sec¬ 
ondly, CaBER and other techniques based on capillary thin¬ 
ning of filaments are well established for viscoelastic fluids 
such as polymer solutions where liquid bridges thin exponen¬ 
tially in time, creating long-lived, slender and almost cylin¬ 
drical filaments.® This permits the use of a simple stress- 
balance equation to extract the viscoelastic fluid stress at the 
necking plane from just a measurement of the radius i? as a 
function of time t. For fluids with little or no elasticity, di¬ 
rectly calculating the viscosity through the stress-balance has 
thus far been shown to be practicable again only for very vis¬ 
cous fluids where a cylindrical filament can form towards the 
final stages of breakup.® Thin cylindrical filaments do form 
close to break-up, but imaging these require the combination 
of ultra-fast and very high resolution photography.® 

The first of these problems was overcome in a method de¬ 
veloped by Bhattacharjee et < 2 /.wherein the liquid-bridge 
is created and stabilized against capillary forces initially by 
power input from surface acoustic radiation (Fig. [^. In our 
experiments, a 20 MHz waveform generator (33220A, Ag¬ 
ilent) was used to generate a surface acoustic wave (SAW) 
burst triggered by a second signal generator (WF1966, NF 
Corporation). The latter delivered a sinusoidal signal to an 
RF power amplifier (411FA, ENI), providing a fixed fre¬ 
quency and amplitude signal near the 36.7 MHz resonance 
frequency of the SAW substrate. An arrangement of curved 
inter-digitated transducers (IDTs) bonded to a piezoelectric 
substrate® focussed Rayleigh waves towards the point where 
a sessile droplet (1 /il; approximately 1 mm dia.) rests (Fig.[^ 
(b)). Energy from the SAW leaks into the droplet causing re- 



Fig. 1 (a) Schematic of experimental setup (b) Curved 
inter-digitated transducers for focussing SAW into a sessile droplet 
(c) Time-lapse images of the formation of a liquid bridge after 
ejection of a jet due to actuation of a sessile drop by SAW; the 
red-line indicates duration of the SAW pulse. 


circulation and bulk motion, leading to an elongated liquid 
that bridges a gap to an opposing parallel surface located 1.5 
mm away from the SAW substrate. The opposing surface was 
coated with teflon and was partially fouled to ensure that the 
jet adheres to the surface but does not spread. The SAW burst 
was ended after 1.5 ms which was found to be sufficient to 
create stable liquid bridges in all our samples. The liquid 
bridge then thins under the action of capillary forces, gener¬ 
ating an extensional flow at the necking plane (Fig.[2(c)). The 
motion of the entire liquid bridge was captured using a high¬ 
speed camera (Photron SA5; 62000 fps; image size: 1.35 mm 
X 2.14 mm (192 x 304 pixels)) with a long-distance video mi¬ 
croscope attachment (K2/SC, Infinity). The set-up is lit by a 
single FED lamp placed behind the filament. The radius of the 
neck in each image frame was obtained using standard image- 
analysis techniques. Initial transients were discarded in each 
case until the neck attained a diameter of 50 pixels (0.352 ± 
0.007 mm); this was taken as the initial time (t = 0) for all 
samples. 

We avoid problems with imaging thin cylindrical filaments 
close to break-up altogether by using neck radius data during 
the early stages of the thinning while the axial curvature is still 
large. This is furthermore advantageous when handling parti¬ 
cle suspensions since the effect of particle interactions with 
the air-liquid interface on the overall dynamics at the neck 
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is weaker when the neck radius is large compared to parti¬ 
cle dimensions. Rescaling the governing equations for slender 
but non-cylindrical Newtonian liquid-bridges^ with the mid¬ 
filament radius Rq ait = 0 (as defined above) and the Rayleigh 
time-scale, 





( 1 ) 


the decay of the rescaled neck radius R"" = R/Ro with rescaled 
time = t/rR is parameterized by the dimensionless volume 
and aspect ratio of the liquid bridge, and the Ohnesorge and 
Bond numbers 
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( 2 ) 


of the regression curve through the data in Fig. is inspired 
by the analytical solution to the inertio- viscous c apillary stress 
balance for a cylindrical fluid filament 


^pR^ = {2X-l)^ + ^^, (3) 

where X is the ratio of the instantaneous axial tension in the 
liquid bridge and 27rjR. If X is a constant, the algebro- 
differential equation above can be treated as a quadratic equa¬ 
tion above for R at any R and t. The negative root of the 
rescaled equation. 


dR^ 

dR 


6 Oh 


(Vl + CR* 



(4) 


where, p, 7 and 7 are the density, surface tension coefficient 
and shear viscosity respectively of a fiuid sample and g is the 
gravitational acceleration. In our experiments, droplet vol¬ 
umes, Rq and the separation between bridge surfaces were 
kept constant. The Ohnesorge number for our samples varied 
by almost three orders of magnitude, whereas the variation in 
Bo was relatively much smaller (around 33 %). Therefore, the 
Bond number can be assumed to be relatively relatively con¬ 
stant. Under such conditions, the rescaled time taken for 
the filament to neck to half its initial radius (i.e. the time to 
R"" = 1 / 2 ) is predominantly governed by Oh. 



Oh 

Fig. 2 Variation of liquid-bridge half-times with Oh: the curve is a 
linear least-squares fit of Eqn. ([^ through the data. 

Figure shows half-time data obtained for a range of 
glycerol-water and sucrose-water mixtures whose properties 
(Ps, 7 and p) were independently measured. Surface tension 
of suspensions were measured using a tensiometer (Analite; 
Selby Scientific). Shear viscosities were measured with a 
Haake Mars (Thermo Scientific) shear rheometer. The form 


where C = (16X - 8)/144 can be integrated with R"" = 1 sls 
the initial condition to obtain 



Strictly speaking, X itself depends on time; it has been 
shown however that, if the neck diameter is significantly 
smaller than the two drops at the end plates, the axial filament 
profile and its dynamics can be well approximated as being 
self-similar. As mentioned earlier, we begin taking measure¬ 
ments of R only after the neck is well formed with a radius 
about half the drop radius at the end-plate. Under such condi¬ 
tions, a constant value of X has been shown to lead to accu¬ 
rate predictions. The values of X given by similarity solutions 
under different conditions have been summarized by McKin¬ 
ley and TripathylEl. When inertia is important and Oh « 1, 
X = 0.5 912 is t he value most likely to be observed in an ex- 
perimentlEEIIIil and hence 0.7135 as Oh ^ 0. When 

viscosity is dominant on the other hand, X = 0.7127^3121 and 
^ 1/2 ^ 7.0522 Oh as Oh ^ 00 . We therefore fit the experi¬ 
mental data obtained for intermediate Oh values with the fol¬ 
lowing rational function: 

* 0.7135+ ii: Oh + 7.0522 Oh^ 

‘7 =-TToh-' ® 

and with K = 14.7 ± 0.2 obtained by linear regression. The 
inverse of the function above is then used as a calibration curve 
to convert observed values of into Oh and then further into 
fj = 3rj with the definition of the Ohnesorge number. The good 
agreement with experimental data of this curve derived from 
the analytical solution indicates that the effect of any residual 
flows generated by the initial SAW irradiation is negligible. 
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3 Modeling 

3.1 Rheology 

The concentration dependence of viscosity of passive sus¬ 
pensions below the isotropic-nematic transition is often mod¬ 
eled on the basis of the mean-field assumption that when a 
new particle is added into the free-volume of a suspension, 
it “sees” the surrounding suspension as a homogeneous fluid 
of higher viscosity than the original suspending medium®^. 
The incremental infiuence of this particle on the overall vis¬ 
cosity is assumed to be independent of the suspension den¬ 
sity itself, and is therefore, in extensional fiows, equal to the 
intrinsic viscosity coefficient of a dilute suspension, [fj] = 
lim^^o “ 3 77s)/ (3 ? 7 s 0), where is the extensional viscos¬ 
ity of the suspension, r]s is the shear viscosity of the suspend¬ 
ing medium, and (j) = nvp is the volume fraction of particles, 
each of volume Vp and at a number density n. T his ar gument 
leads to the Krieger-Dougherty equation for the 

relative extensional viscosity of a non-dilute suspension. 



where <prn is the maximum volume fraction at which the 
steady-state viscosity diverges and beyond which suspensions 
behave as slurries with significant solid-like characteristics. 

Recent studies by Haines et a/.^and Saintillanl^have mod¬ 
eled the statistics of an axisymmetric self-propelled particle 
that stochastically changes its swimming direction while mov¬ 
ing in an externally imposed homogeneous extensional fiow, 
and then obtained an expression relating [fj] to motility char¬ 
acteristics and strain-rate. We summarize here Saintillan’s 
analysis for a slender rod to explain the physical significance 
of the key parameters in the expression for [rj] in an active sus¬ 
pension. The particle contribution to the fiuid stress tensor is 
proportional to the orientational average of its total hydrody¬ 
namic dipole which is the sum of dipoles arising from propul¬ 
sion, the externally-imposed fiow, and thermal fiuctuations. 
Besides translating with a mean speed U, living swimmers 
also autonomously change direction in an apparently random 
manner which is modeled through a non-Brownian rotational 
diffusivity, Dj.. Passive particles (e.g. dead swimmers) on the 
other hand only have a diffusivity o that is solely due to 
thermal fiuctuations. The suspending medium in turn expe¬ 
riences an equal and opposite reaction from the particle. In 
a dilute suspension undergoing a homogenous fiow, the parti¬ 
cle contribution Xp to the fiuid stress tensor thus arises from 
an average of hydrodynamic dipole moments induced at the 
particle location by these forces i.e. 

Tp = n [(S^> + {S^> + (S^’)] , (8) 

where angular brackets represent an ensemble average over 
the distribution of particle orientations p and n is the particle 


number density. In the Stokesian regime, the fiow-induced 
dipole is proportional to the rate-of-strain tensor E, with a 
shape-dependent proportionality constant A. The propulsive 
dipole strength is represented by another shape-dependent 
constant a. The particle stress is thus derived to be 

Xp = nA ^{pppp>- i(pp)6j:E 

+ (3nfcBT+ncr) ^{pp>-i6j , (9) 

where E and 6 are the rate-of-strain and identity tensors, re¬ 
spectively, and /cb is the Boltzmann constant and T is the 
absolute thermodynamic temperature of the suspension. The 
orientational distribution of a slender rod is governed by a bal¬ 
ance between frictional interactions with the medium that ro¬ 
tate the particle and the stochastic kicks that tend to make the 
distribution isotropic. The resulting Fokker-Planck equation 
has a closed form analytical solution for an extensional fiow at 
steady-state. The averages in the equation above can hence be 
obtained for a uniaxial extensional fiow of strain-rate e. Defin¬ 
ing the particle contribution to the suspension extensional vis¬ 
cosity 77p = {xp^zz - ^p,rr)I where the principal stretching 
direction is along the ^-axis while r is the transverse radial 
coordinate, the specific extensional viscosity ratio is given by 


^Vs 


nA 

‘^Vs 



3 Dr I /cbT a 
e \ ADr 3ADr 



( 10 ) 


where /3 is a shape-dependent constant that arises in the cou¬ 
pling of the rotational rate of the particle with the straining 
motion of the fiuid around it.^ The function M depends on 
P e and is discussed shortly. We first identify the key dimen¬ 
sionless model parameters and relate them to parameters that 
are physically measurable. 

For slender rods of hydrodynamic length and diameter 
d it is known that^S 


^Vs Ll 
6 1n(2Ljd) ■ 


( 11 ) 


In order to apply the analysis above to suspensions of mi¬ 
crobes, we must assume that the frictional characteristics of 
more complex axisymmetric shapes can be represented by an 
effective cylinder. The total length L for a fiagellar swimmer 
is the sum of the lengths of the head and the tail when com¬ 
pletely stretched out. A living swimmer does not of course 
take a fully-stretched shape, and its effective hydrodynamic 
length Lh while swimming will in general be less than L. 
Defining the ratio A = L^jL, the dimensionless equivalent of 
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A is obtained as 


A _ ttA^ 

2r]sVp 12 Vp \n{2XL/d) 


( 12 ) 


The propulsive thrust generated by the inertialess swimmer 
and the corresponding frictional resistance to swimming are 
equal and opposite to one another, but act at different loca¬ 
tions on the axis. A scaling argument for the strength of this 
propulsive dipole exerted by particle swimming at a speed U 
suggests (7 ~ r]sL^U. Defining a geometry-dependent pre¬ 
factor a such that a = {t: l2) arjs L‘^ U, we obtain a dimen¬ 
sionless equivalent of a as 


a a\n(2\Lld)U 

^ ~ 3ADr X^LDr ■ 

The sign is positive for pullers and negative for pushers, while 
for passive particles, ao = 0. Brownian torques due to thermal 
fiuctuations tend to relax the orientational distribution of par¬ 
ticles to its isotropic equilibrium state; their relative strength 
in an active suspension is represented by 


__ /cbT _ 6/cbT In (2 AL/d) 

^ A Dr 7rr]s Dr 

The Fluctuation-Dissipation Theorem (FDT) is valid for pas¬ 
sive particles, in which case the rotational diffusivity for slen¬ 
der rods is 


^ kBT 3teTln(2(AL)/d) 

^r,o = =-TVTAs-• 

2 A 7rr]s{XLy 

Therefore, 70 = 2 for dead cells; no such simple relation is 
currently available for 7 in active systems. 

The Peclet number — Pe = ejDr (or e/D^^o)— quantifies 
the relative ability of the extensional flow of strain-rate e to 
reorient a particle against random switches in swimming di¬ 
rection. Using the definitions of the dimensionless parameters 
in Eq. ( p^ , the intrinsic extensional viscosity is obtained as 


[7] = lim 


% 




= X 




3 , 

( 1 \ 

/ 1 \ 

+ — 


M- - 

Pe ' 

1 /?/ 

V 3 / 


(16) 


The function M(x) is weakly monotonic in x and has the form 

1 1 


M{x) = 


(17) 


2xD{x) 2x‘^ ’ 

where D{x) = exp(-x^) exp(^^) dy is Dawson’s integral. 

In Eqn. above, x = y^3/3Pe/4. As Pe ^ 0 therefore, 
M = 1/3 + 4/45x2 + O(x^) = 1/3 + pFe/lb + O(Pe^), and 


when Pe ^ 00 , M{x) = 1 - x ^ + 0{x ^) = 1 - 4/(3/5 Pe) + 
0(Pe“^). 

The model above for [fj] requires, besides the solvent con¬ 
ditions kBT and r]s^ independent measurement of the follow¬ 
ing parameters: the average geometric characteristics of the 
swimmers, L, d, Vp\ their motility characteristics U and Dr 
(or Dr^ 0 for dead cells); and the hydrodynamic ratios A and a. 
The constant (3 is taken to be unity, whic h is appropriate for 
slender particles of large aspect ratioOut of these, the sol¬ 
vent parameters are obtained by standard techniques and the 
particle size characteristics L, d and Vp are obtained through 
microscopy. As noted above, the value of Dr^ 0 is effectively 
set by invoking the FDT and using 70 = 2. We discuss be¬ 
low the estimation of the motility parameters, U and Dr. It is 
more difficult to directly determine the ratios A and a, which 
are therefore treated as free parameters to be obtained by com¬ 
paring model predictions to experimental data for Keep¬ 
ing in mind their physical significance, we assess later if the 
values for A and a obtained thus are plausible. 


3.2 Active diffusivity 

In the absence of fiow, the trajectory of an active particle 
swimming in with a speed U along its principal axis with ran¬ 
dom reorientations of that axis can be modeled by a pair of 
Langevin equations for its instantaneous position r(t) and ori¬ 
entation. The rescaled mean-squared displacement for such a 
particle in two dimensions is^ 

= 4 ^ 1+2 (i-l + e-^) , (18) 

where i = Dr t and ^ = Dt DrlU^. The rescaled translational 
diffusivity ^ controls the shape of the MSD-versus-t curves: 


MSD 


I (4^ +t)t, iff « 1 
[(4^ + 2 )f, iff»l. 


(19) 


If translational diffusion is dominant and ^ » 1, then MSD 
4^t always. When 0 < ^ < 1 on the other hand, one ob¬ 
serves diffusive-ballistic-diffusive behaviour. At short times 
when t « 4^, MSD 4^t. At long time scales too, the be¬ 
haviour is diffusive and MSD (4^+2) t. Therefore, if the ef¬ 
fective long-time translational diffusivity Deff is defined such 
that MSD = 4Defft, then Dq^ = Dt + _{U^l2 Dr). When 
4^ « t < 1, the motion is ballistic with MSD P. Values of 
Dr, Dt and U can be obtained from shift-factors by shifting 
experimental MSD-versus-t data to match one of the dimen¬ 
sionless theoretical curves on a log-log plot. 
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4 Results and discussion 

Experimental MSD data was obtained using suspensions of 
very low volume fractions with cell number densities smaller 
than L~^. We find that this data for all the three organisms 
resembled the ^ = 0 curve, indicating that translational dif- 
fusivity during swimming has a negligible influence (Fig.|^. 
The values of swimming speed and diffusivity thus extracted 
are reported in Table The value of for algae is com¬ 
parable to those observed by Rafai et. in another blue- 

green algal species C reinhardtii: they measured an average 
speed ofU = 40 /im/s, and a long-time effective diffusiv¬ 
ity of ^ 995 /im^/s, which corresponds to a rotational 
diffusivity of Dr ~ 0.8 s“^ assuming that the true transla¬ 
tional diffusivity is negligible in that species as well. For E. 
coli, analysis of the 3D cell-tracking data suggests Dr ~ 3.5 
s-i[34|^ although a much lower value of Dr = 0.057 s“^ was 
reported by Drescher et al^ Some of this variability can be 
attributed to differences in strains and the media and protocols 
used for cell-tracking. In addition, there are also significant 
differences in the analyses of the run-and-tumble motion of E. 
coli cells to assign effective diffusivities. Although computer- 
aided analysis of sperm motility is well established in mam¬ 
malian reproductive biology,®! to our best knowledge, mean- 
square displacements have not been analyzed to determine ef¬ 
fective diffusivities. The values of U extracted by fitting the 
theoretical prediction in Eqn. ( p^ through the mean-squared 
displacement data for the algae and bacteria are close to the 
values measured by direct observation (45 and 4.7 /xm/s, re¬ 
spectively). In the case of sperm however, the extracted value 
is significantly different from the direct measurement of an av¬ 
erage speed 70 yum/s for the motion of the sperm head which is 
comparable to average hea d spe eds of around 100 - 150 /im/s 
in hyperactive mice sperm.®^. This difference could be due 
to the fact that the sperm head oscillates about the mean tra¬ 
jectory of the cell.®! 


Table 1 Geometric and motility characteristics 



Algae 

Bacteria 

Sperm 

L (fim) 

36 

13 

95 

d (/xm) 

2.8 

1 

3 

Ljd (/xm) 

12.9 

13 

32 

U (/xm/s) 

30 

5.5 

31 

Dr (1/s) 

3 

10 

0.8 


Figure [^presents the key results in our study. The top panel 
in the Figure shows the evolution of the neck radius during 
capillary thinning for the highest volume fractions studied in 
each species; data at lower volume fractions is presented in 
the Supplementary Information.^ We observe that capillary 



Fig. 3 Comparison of experimental data for the mean-squared 
displacement for algae, E. coli and mouse sperm with predictions 
for various values of the relative translational diffusivity ^; 

MSD = MSD {DrjUf and i = Dr t 


thinning of samples of live algal suspensions progresses more 
slowly than of samples of dead cell suspensions at the same 
volume fraction (top-panel, Fig.|^. The opposite behaviour is 
observed in bacterial and sperm suspensions. This is in line 
with the expectation that motility of pullers like D. tertiolecta 
tends to increase the viscosity whereas in pushers such as E. 
coli and sperm, viscosity decreases. The middle and bottom 
panels in Fig/|^ shows that the differences between the radius- 
versus-time data, and thus between the viscosities, of live and 
dead cell suspensions. These differences in E. coli are small 
but significant relative to experimental uncertainty. In com¬ 
parison, the differences are larger in algae over a similar range 
of volume fractions. Data for sperm could only be obtained at 
low concentrations due to the small sample volumes collected 
from mice and limitations in concentrating them further. How¬ 
ever, even at these low concentrations, the effect of motility 
on viscosity appears clear.The difference in viscosity between 
live and dead cell samples appears to increase with volume 
fraction for all three species. 

To use the model to understand these trends, we first need 
to determine the strain rates in the experiments. As is well 
known, the instantaneous strain-rate e = -2dhiRldt at the 
necking plane in CaBER experiments cannot be directly con¬ 
trolled but are determined by the liquid bridge dimensions and 
fluid properties. The radial decay in all our samples was ob¬ 
served to be approximately linear when t < ti/ 2 . The in¬ 
stantaneous strain-rate in this observation period for linear 
radial decay increases from l/ti /2 to 2 /ti/ 2 . We can thus 
estimate the average strain-rate from the obesrved ti /2 as 
^avg ^ ^ 1/2 /ti /2 = 1 . 4 /ti /2 since the the Hencky strain at ti /2 
for a linear radial decay can be shown to be e 1/2 = 2 In 2. It is 
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Fig. 4 Top panel: radial decay observed during capillary thinning for the most concentrated algal, bacterial and sperm suspensions studied; 
middle panel: variation of dimensionless capillary thinning half-times with particle volume fraction; bottom panel: comparison of predictions 
of relative extensional viscosity = fy/(3?7s) with experimental values extracted from measured = ti/ 2 l^R- In B and C, up- and 

down-triangles represent data for live and dead cell suspensions, respectively. Continuous curves in the bottom panel have been obtained with 
values of a and A as shown on the plots whereas dashed curves are predictions with both these parameters set to unity. 
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therefore clear from the systematic increase in ti /2 with con¬ 
centration (Fig. |4]B (i — iii)) that the average strain-rate Cavg 
decreases with concentration as the increasing viscosity slows 
down capillary thinning. The strain-rate at any concentration 
in the range covered in the experiments was obtained using a 
cubic polynomial fit through the experimental Cavg-vs-^ data.^ qj 
Values of Cavg in our experiments ranged from about 250 s“^ 
at higher concentrations to 2000 s“^ for pure buffers. 


Table 2 Model parameters 



Algae 

Bacteria 

Sperm 


Free parameters 


a 

15 

1 

30 

A 

0.6 

0.3 

0.7 

pm 

0.7 

0.7 

0.7 


Calculated parameters 


X 

5.2 

6.7 

100 

7 

4.5 X 10“^ 

1.2 X 10"^ 

5.4 X 10“^ 

a 

40 

-2.7 

-97 

F(pN) 

46 

0.6 

260 

Dr, 0 (1/s) 

6.8 X 10“^ 

6.1 X 10“^ 

2.2 X 10“^ 


Values of model parameters and other derived quantities for 
the three species are listed in TableThe ratio A of the hydro- 
dynamic length to the total end-to-end length L is a parameter 
that is estimated by visually bringing the predictions for dead 
cell suspensions into agreement with the experimental data. 
Since we are interested in order of magnitude estimates for this 
and the other free parameters a and (pm in the model, a more 
rigorous least-squares fit was not pursued. The estimation of 
the parameter A is not independent of the maximum packing 
fraction <pm- While <pm ^ 0.63 is often used for spheres, there 
is no consensus on its value for anisotropic rod-like particles. 
It has been shown that pm for such systems depends on par¬ 
ticle aspect ratio for dilute systems.® More concentrated sus¬ 
pensions of rod-like particles undergo an isotropic-to-nematic 
transition w hich is expected to occur ai p ~ 0 {{Lld)~^) at 
equilibrium.®^ The aspect ratios (based on the total end-to- 
end length; Table of all three species are quite large. The 
concentrations of the algal and bacterial suspensions in our ex¬ 
periments are thus possibly well above the isotropic-nematic 
transitions for those systems. Although the volume fractions 
of the sperm suspensions are low, significant orientational ef¬ 
fects due to interparticle interactions cannot be ruled out due 
to their large aspect ratios. A value of pm = 0.7 gives reason¬ 
able agreement of model prediction with dead cell data for all 
three species. With A and pm determined in this manner, a is 
estimated to match predictions for live cells with experimental 
data for 

The predictions of in Fig. were generated using the 
KDE after first calculating the intrinsic viscosities [rj] with 



0 


Fig. 5 Experimentally observed variation of Pe with volume 
fraction in live samples; the curves have been obtained by 
polynomial regression of strain-rate data.^ 


Eqn. ( p^ for the ranges of Peclet numbers observed for live 
and dead cells. The Peclet numbers calculated as Pe = Cavg Dr 
with the rotational diffusivities for live cells in Table [T] are 
plotted in Fig.|^ The Peclet number governs the competition 
between fiow and diffusion and the values of Pe » 1 in Fig.[^ 
indicate that the fiow completely overcomes the tendency for 
diffusion to make the orientational distribution isotropic and 
that particles are aligned along the principal stretching axis of 
the unaxial extensional fiow. In such an aligned state, the rel¬ 
ative contributions of the flow-induced and propulsive dipoles 
to the total stress and hence the viscosity is determined by the 
magnitude of the parameter a in an active suspension. In fact, 
in the case of dead cells for which this parameter is zero, the 
fiow wins by default. Moreover, the much smaller rotational 
diffusivities (Table of dead cells result in even larger Peclet 
numbers than in live cell suspensions at comparable strain 
rates. The values of [fj] for dead cells in each species are thus 
nearly identical to the limiting value [7j]°°2x/3 contributed 
solely by the flow-induced dipole as Pe ^ oo. In active sus¬ 
pensions on the other hand, even when cells are highly aligned 
by the fiow when Pe » 1, the contribution of propulsive forces 
can be relatively large and [rj] is significantly different from 
[rj] °°. This is the principal reason behind the significant dif¬ 
ferences observed between live and dead cell samples in Fig.[^ 
despite the large strain rates and Pe. Moreover, with increas¬ 
ing cell volume fraction, not only are there more cells per unit 
volume contributing propulsive stresses, the strain-rate and Pe 
numbers decrease (Fig.|^ boosting the relative contribution of 
active dipole. 

The effect of activity on Fig. [^appears strongest in 

algae because of the combination of a large propulsive dipole 
and particle concentrations. The dipole is much weaker in E. 
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Fig. 6 Effect of flow strength on intrinsic viscosity of live cell 
suspensions of sperm (continuous), algae (dashed) and bacteria 
(dot-dashed curve): the dark portions of the curves represent the 
experimentally observed ranges of Pe values; the symbols represent 
the zero-Pe limit, /[^]°° obtained with the model parameters in 
Table[^ 

coli pushers, leading to a smaller effect of activity, although 
they have the largest swimming diffusivity and lowest Pe val¬ 
ues among the three species. In contrast, a in sperm appears 
to be so large that their live suspensions show a clearly mea¬ 
surable effect of motility even at very low cell densities and 
despite very large Pe values. Interestingly, the negative values 
of a for the pushers E. coli and sperm are predicted to lead 
to negative [f]] at sufficiently low strain-rates (Fig. [^, which 
will result in live cell suspensions being less viscous than the 
suspending medium. Although such low strain-rates cannot 
be practically realized with extensional rheometers currently, 
they are accessible in shear rheometers. This intriguing ef¬ 
fect has been demonstrated by Gachelin et al^ with E. coli in 
a microfiuidic device. This effect is likely to be much more 
dramatic with sperm suspensions: the parameters we have ob¬ 
tained for sperm lead to a very large negative zero-Pe limit 
[v]^ = X (2/9 + P/lb (7 a - 1/P )) = -1.9 x 10^. 

How reasonable are the values (in Table of the key free 
parameters: the prefactor a in Eqn. ( p^ for the active dipole, 
and the ratio A of the hydrodynamic length to the total head- 
tail length of a swimmers? If both these parameters are set 
to unity, the values of the activity a are too low for all three 
species and predictions for live and dead cell suspensions are 
virtually indistinguishable, and these predictions are substan¬ 
tially different from the experimental data (dashed curves in 
Fig.gca — iii)). The values of A < 0(1) required to obtain 
agreement with experiment appear reasonable given that flag¬ 
ella are never fully stretched out in swimming cells. Hydrody¬ 
namic interactions between the head and tail cause the param¬ 
eter a to be strongly geometry dependent. This can be seen by 


modeling a fiagellar swimmer modeled as an asymmetric rigid 
dumbbell with Stokeslets of different hydrodynamic radii lo¬ 
cated at the head and tail respectively (Appendix A).l^We ob¬ 
tain for such a dumbbell a = 12 {Lht/L)‘^ / [{Lht/dh) - 3/2], 
where L^t is the distance between the head and tail centres 
and Oh is the hydrodynamic radius of the head. Although 
finite-size corrections can be expected to modify the singu¬ 
lar behaviour when Lht/dh = 3/2, this result suggests that a 
can vary over a wide range depending on the ratio L^t and 
Viewed in this light, the variation in values of a for the three 
very differently shaped cells appears plausible. 

It if further possible to estimate propulsive forces (using 
the definition of a in Eqn. and its values in Table 
as F = cr/Lh. Bayly et alP^ analyzed the fiagellar stroke 
of a single C. reinhardtii using resistive force theory to esti¬ 
mate an average power dissipation of about 5 fW, which when 
combined with an average speed of around 40_um/s in that 
species,® yields F = 125 pN. Drescher et a/.® estimated a 
propulsive force of F = 0.43 pN directly in a wild-type F. 
coli cell by measu ring the flow field around it during a straight 
run. Schmitz et < 2 /.® used a microprobe to measured the force 
required to stall the motion of a bull sperm fiagellum to be 
F = 250 pN. These values from direct single cell or fiagellum 
measurements are comparable to the values of F in Table [^. 

5 Conclusions 

Our experimental observations with a surface-acoustic-wave 
driven microfiuidic rheometer indicate that particle motility 
has a clearly measurable infiuence on the rheology of suspen¬ 
sions. Capillary thinning of liquid bridges proceeded more 
slowly in suspensions of algal pullers than those of dead cells 
at the same volume fraction, whereas bacterial and sperm 
pushers tended to hasten thinning. The difference in the ef¬ 
fective viscosity between suspensions of live and dead cells 
was found to systematically increase with concentration. Pre¬ 
dictions with a model that combined the Krieger-Dougherty 
equation for the relative viscosity of suspensions with an equa¬ 
tion derived by Saintillan^ for the intrinsic extensional vis¬ 
cosity of active suspensions were found to be in good quali¬ 
tative agreement with the experimental observations. Our re¬ 
sults show that the propulsive dipole even in weak swimmers 
such as F. coli can contribute significantly to fiuid stresses 
even at high Peclet numbers when flow dominates over swim¬ 
ming noise. While parameters such as the properties of the 
suspending medium, the physical dimensions of the cells, and 
their swimming speeds and effective rotational diffusivities 
could be determined from independent measurements, model 
predictions depend crucially on hydrodynamic details such as 
the effective hydrodynamic length of the swimmer and the 
contribution of hydrodynamic interactions between the head 
and tail to the propulsive dipole moment. Values of these pa- 
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rameters required to match predictions with experimental data 
lead to estimates for propulsive forces in the three species that 
are in line with data in literature on direct force measurements 
with single cells or flagella. This suggests that it is possi¬ 
ble to develop a better models for the nonlinear rheology of 
microswimmer suspensions by combining more accurate al¬ 
ternatives for the Krieger-Dougherty formulation describing 
concentration effects with further refinements of the dilute so¬ 
lution theoryl^E] for the viscosity of active suspensions that 
more accurately reflect the hydrodynamic characteristics of 
swimming particles and their propulsion. Such models open 
up the possibility of routinely and precisely extracting propul¬ 
sive forces from rheometry of active suspensions. 

6 Acknowledgments 


to the translational motion: 

F, = -Ct(u + n.F;,) +F^, (20) 

where (t = 6 tt is the drag coefficient of the tail and O 
is the tensor describing the hydrodynamic interaction tensor 
between the head and tail such that the velocity perturbation 
at the tail due to the hydrodynamic force on the head is 
-O • Fh^. The head only experiences frictional resistance, 
and 

Fh = -Ch(u + Fl^Ft). (21) 

The total hydrodynamic force on the inertialess swimmer must 
be zero in the absence of external forces: 

Fh + Ft = 0. (22) 
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A Self-propelled asymmetric dumbbells 


The three equations above can be solved for F^, Ft and Fp, 
given the frictional characteristics ((h, Ct the ve¬ 

locity u. Assuming that all forces and the velocity are directed 
along the axis of the dumbbell, if Fh, Ft and Fp, and u are the 
axial force and velocity components, we obtain, 

Ft = -Fh = —(23) 

(24) 



where U = 1/(4 tt L^) for Oseen-Burgers hydrodynamic in¬ 
teraction, and Lp is the distance between the head and tail cen¬ 
tres. The magnitude of propulsive dipole is then 


kl = \Fh\Lp = 


12 {n 12) rjsLl !«! 
{Lpl ah) - (312) ■ 


(25) 


With a defined as before i.e. |(t| = (TTl2)aris |m| where L 
is the total end-to-end length of the swimmer, we obtain 


l2{Lp/L)^ 

{Lp/ah) - (3/2) • 


(26) 
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Fig. 7 Self-propelled asymmetric dumbbell model of a 
microswimmer with a flagellar tail 

A swimmer is modeled as a rigid asymmetric dumbbell with 
Stokeslets of hydrodynamic radii ah and at located at the hy¬ 
drodynamic centres of the head and tail respectively!^. The 
rigid swimmer translates with a velocity u relative to the un¬ 
perturbed ambient fluid as a result of the actuation of the tail. 
The tail generates a thrust of F^. The total hydrodynamic 
force on the tail is the sum of Fp and its frictional resistance 
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